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We study the UV behaviour of actions including integer powers of scalar curvature and even powers 
of scalar fields with Functional Renormalization Group techniques. We find UV fixed points where 
the gravitational couplings have non-trivial values while the matter ones are Gaussian. We prove 
several properties of the linearized flow at such a fixed point in arbitrary dimensions in the one-loop 
approximation and find recursive relations among the critical exponents. We illustrate these results 
in explicit calculations in d = 4 for actions including up to four powers of scalar curvature and two 
powers of the scalar field. In this setting we notice that the same recursive properties among the 
critical exponents, which were proven at one-loop order, still hold, in such a way that the UV critical 
surface is found to be five dimensional. We then search for the same type of fixed point in a scalar 
theory with minimal coupling to gravity in d = 4 including up to eight powers of scalar curvature. 
Assuming that the recursive properties of the critical exponents still hold, one would conclude that 
the UV critical surface of these theories is five dimensional. 
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I. INTRODUCTION 

In [l[ scalar-tensor theories were studied where the purely gravitational part was given by the Einstein-Hilbcrt 
action. Here we generalize those results by including higher curvature terms. The main aim of our analysis is to 
understand if gravity remains asymptotically safe 0, HJ under the inclusion of some matter component. Results 
about the renormalizability of gravity can depend crucially on the inclusion of matter. Already in the first one-loop 
calculations 0, [E[ it was shown that pure gravity is one-loop renormalizable but becomes one-loop nonrenormalizable 
in the presence of matter. In the context of the search for asymptotic safety, it was shown in Q that the position and 
even the existence of a nontrivial gravitational fixed point in the Einstein-Hilbcrt truncation is affected by the presence 
of minimally coupled matter fields. In Q we showed that this is also true when higher derivative gravitational terms 
are present. In [l|, [&, 9] the effect of gravity on scalar interactions was studied, assuming the Einstein-Hilbcrt action 
for the gravitational field. In [l], Q it was shown that a nontrivial fixed point exists, where the purely gravitational 
couplings arc finite while those involving the scalar field vanish. This is called a "Gaussian matter fixed point" 
(GMFP). In the present paper we extend these results by considering an interacting scalar field coupled to a class 
of higher derivative gravity theories which had been studied previously in @, [l(| ■ We ask whether the scalar matter 
contribution is able to alter the results of the purely gravitational part considerably. In Q the addition of minimally 
coupled matter components to i? 2 -gravity (including all possible curvature invariants up to quadratic order) showed 
that the nontrivial fixed point structure is maintained in that case. We will see that this is largely the case here too, 
but since we now consider interacting scalars we will find that the dimension of the critical surface increases. 

There are clearly many possible applications in cosmology. Early work in this direction has been done in (ll| , using 
the beta functions of pure gravity. Taking scalar fields into account could have significant renormalization group 
running effects in inflation. Without the necessity of asymptotic safety, in effective field theory calculations the beta 
functions derived here could be useful e.g. for inflation [lq . or in models where the Higgs field is used as the inflaton 
field Applications in the IR are possible, for example along the lines of [l4| or the much discussed modified 

theories of gravity with some action based on different functional forms of the Ricci scalar (see e. g. (Hi). We 
mention that the appearance of a scalar field in the low energy description of gravity has also been stressed in [r|. 
For a FRGE-based approach to that issue see also [l7| • 

As in [l| the analysis is based on a type of Wilsonian action Tk called the "effective average action" depending on 
an external energy scale k which can be formally defined by introducing an IR suppression in the functional integral 
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for the modes with momenta lower than k. This amounts to modifying the propagator, leaving the interactions 
untouched. Then one can obtain a Functional Renormalization Group Equation (FRGE) [18| for the dependence of 
T k on k, 



d t T k = -STr 
2 
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where t = log(/c / ko) . $ are all the fields present in the theory. STr is a generalized functional trace including a minus 
sign for fcrmionic variables and a factor 2 for complex variables. TZ k is the regulator that suppresses the contribution 
to the trace of fluctuations with momenta below k. As the effective average action contains information about all the 
couplings in the theory, the FRGE contains all the beta functions of the theory In certain approximations one can 
use this equation to reproduce the one-loop beta functions, but in principle the information one can extract from it 
is nonperturbative, in the sense that it does not depend on the couplings being small. 

A Quantum Field Theory (QFT) is asymptotically safe if there exists a finite dimensional space of action functionals 
(called the ultraviolet critical surface) which in the continuum limit are attracted towards a Fixed Point (FP) of the 
Renormalization Group (RG) flow. For example, a free theory has vanishing beta functions, so it has a FP called the 
Gaussian FP. Perturbation theory describes a neighbourhood of this point. In a perturbativcly renormalizable and 
asymptotically free QFT such as QCD, the UV critical surface is parameterized by the couplings that have positive or 
zero mass dimension. Such couplings are called "renormalizable" or "relevant" . Asymptotic safety is a generalization 
of this behaviour outside the perturbative domain. That means that the couplings could become strong. The FRGE 
allows us to carry out calculations also in that regime. 

Whether gravity is indeed asymptotically safe cannot yet be fully answered. However, since the formulation of the 
FRGE by Wetterich [l8[ , many results support this possibility in various approximations [13, [H, H3, HH, H3, [H, [24| , 
for reviews see [25| . Since the first application to gravity the necessary tools have been developed to make the approx- 
imation schemes more reliable including more couplings and studying their UV behaviour. In the approximations so 
far used, gravity has a nontrivial fixed point with a finite dimensional UV critical surface as is consistent with the 
requirements of asymptotic safety. 

The most common approximation method is to expand the average effective action in derivatives and to truncate 
the expansion at some order. In the case of scalar theory the lowest order of this expansion is the local potential 
approximation (LPA), where one retains a standard kinetic term plus a generic potential. In the case of pure gravity, 
the derivative expansion involves operators that arc powers of curvatures and derivatives thereof. This has been 
studied systematically up to terms with four derivatives in (2ll . l22l . |23| . |24| and for a limited class of operators (namely 
powers of the scalar curvature) up to sixteen derivatives of the metric [a, [l^ . In the case of scalar tensor theories of 
gravity, one will have to expand both in derivatives of the metric and of the scalar field. 

In this paper we will study the generalization of the action considered in P, 0, H| and [l(| of the form 



~GF 
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where Sqf is a gauge fixing term to be specified below and S g h is the corresponding ghost action. This action can 
be seen as a generalization of the LPA where also terms with two or more derivatives of the metric are included. 

This paper is organized in the following way. In section [II] we will give the inverse propagators resulting from the 
action ([2]) which have to be inserted into the FRGE to obtain the beta functions. In section IIIII we describe the 
general properties of the GMFP. It is divided into two sub-sections. In section fill Al we show that minimal couplings 
are self-consistent in the sense that when matter couplings are switched off then also their beta functions vanish. 
In section IIII B[ we analyze the linearized RG flow around the GMFP. We find that the stability matrix has block 
diagonal form which allows us to calculate its eigenvalues and eigenvectors in a recursive way. In section IIV Al we 
illustrate the existence of the GMFP and the properties of the RG flow near the GMFP in specific truncations where 
scalar matter fields are coupled nonminimally to gravity, including operators with up to four powers of scalar curvature 
and quadratic in the scalar matter field. In section ITVBI we consider minimally coupled scalar tensor theory including 
operators up to eight powers of scalar curvature and determine the dimensionality of the UV critical surface. We 
conclude in section IVl 
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II. THE FRGE FOR F((f> 2 ,R) 
A. Second variations 

Starting from the action given in eq. ((2]), we expand F((f> 2 ,R) in polynomial form in <fi 2 and R as 

p 

F(0 2 , R) = VbW> 2 ) + V^ 2 ) R + V 2 {4> 2 ) R 2 + W 2 ) R 3 + ■ ■ ■ + V p {^ 2 ) R p = J2 W*) RU 

a=0 



(3) 



In order to evaluate the r.h.s. of eq. ([T]) we calculate the second functional derivatives of the functional given in 
cq. (J2J) . These can be obtained by expanding the action to second order in the quantum fields around classical 
backgrounds g^ — g^ + and <fi = cf> + 6(j>, where </> is constant. The gauge fixing action quadratic in is chosen 
to be 



(4) 



where \v = V M /i M „ — i^V^/i^, G M „ = <? M „ (a + /?□); a, /3, and p are the gauge parameters, we denote □ = V^V^. 

The gauge fixing action eq. flU gives rise to a ghost action consisting of two parts, S g h = S c + 5&. The first part 
S c arises from the usual Fadeev-Popov prodedure leading to the complex ghost fields C M and C^. It is given by 



d d x^C p (a + f3D) 



c„ 



(5) 



The second part Sb arises for (3 ^ and comes from the exponentiation of a nontrivial determinant which requires 
the introduction of real anti-commuting fields which are usually referred to as the third ghost fields [26| , 

Sb=\j dtxyfijbuG^bu. (6) 
These terms are already quadratic in the quantum fields. Then the second variation of eq. @ is given by 



,(2) 



2 / d x V9 



F(tf,R)\-h 2 --h^ 



h^R^ A h pX - h»V^ p h pu + WVV 



| + dF{ f R R) {- hh ^ R n» - \ hUh + 5^ n V + W a h a0 Rl 
bTR^ ■ h al3 R a - 2h^R^ ■ VV ff V 



d 2 F{<j) 2 ,R) 
dR? 



-iWRw ■ ah + v a v p h aP ■ vvv - 2Dh ■ V'VV + nh ■ nh 



d d x^g 



h-<j> 



l - J d d x^g8<j> 



dF{4> 2 ,R) 
d(f) 2 

-□ + 2 



5<f> + 2(j)5<j> 



d 2 F{cp 2 ,R) 
dRdcb 2 



VVhur-nh-h^Rw 



dF^^R) , ~&F{<p,R) 



d<j> 2 



4^ 



d{4> 2 Y 



SqF + Sgh 



(7) 



where □ = V^V^ and h = h)t. Since we will never have to deal with the original metric g^ v and scalar field </>, in 
order to simplify the notation, in the preceding formula and everywhere else from now on we will remove the bars 
from the backgrounds. As explained in [l9j], the functional that obeys the FRGE ([I]) has a separate dependence on 
the background field g^ v and on a "classical field" (gd)^ = g^v + (hd)^, where (h c i) is the Legendre conjugate 
of the sources coupling linearly to {hd)^- The same applies to the scalar field. In this paper, like in most of the 
literature on the subject, we will restrict ourselves to the case when (g c i) uv = g^v and <j> c \ — <j). From now on the 
notation g M „ and <j> will be used to denote equivalently the "classical fields or the background fields. 



B. Decomposition 

In order to simplify the terms and partially diagonalize the kinetic operator, we perform a decomposition of in 
tensor, vector, and scalar parts as in [(j [ic| . 

hpv = h^ v + + V„f M + V^V^cr - ^-g^Ua + ^9^ h ( 8 ) 
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where K^ v is the (spin 2) transverse and traceless part, £ M is the (spin 1) transverse vector component, a and h 
are (spin 0) scalars. This decomposition allows an exact inversion of the second variation under the restriction to 
a spherical background. With that in mind, we work on a d-dimcnsional sphere. For the spin-2 part, the inverse 
propagator is 



ldF((f> 2 ,R) 
2 dR 



□ ■ 



2(d- 2) 
d(d-l)' 



R\- l -F{4> 2 ,R) 



(9) 



where 



5 2 r k 



g^ a g vp ) . For the spin-1 part it is 
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The two spin-0 components of the metric, a and h, mix with i 
3x3 matrix S with the entries 
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resulting in an inverse propagator given by a symmetric 
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As discussed in more detail in [y], to match the trace-spectra of the Laplace-operator acting on h^ u with those 
obtained for the constrained fields after the decomposition, the first eigenmode of the operator trace over the vector 
contribution and the first two eigenmodes of the operator trace over the a contribution have to be omitted. The trace 
over the h and 8<p components should be taken over the whole operator spectrum instead. To handle the mixing of the 
scalar components in an easy way, we subtract first the two first eigenmodes from the complete scalar contribution 
from the matrix S and then add the first two trace modes which should have been retained for h and 8<p. This requires 
to take into account a further scalar matrix B formed by the components of h, <f> and their mixing term. It is given 
by 



B = 



Shh Sh$ 
Sh<j> 5*00 



(12) 



whose trace contribution to the FRGE will be calculated on the first two eigenmodes of the spectrum of the Laplacian. 

Again, in order to diagonalize the kinetic operators occurring in the ghost actions eqs. ((5]) and ([6|), we perform a 
decomposition of the ghost fields C^, and into transverse and longitudinal parts, 



V„C, 6 M = bj l + V M 6, (13) 

After this decomposition, the inverse propagators for the vector and scalar components of the ghost and third ghost 
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C. Contributions by Jacobians 



The decomposition of h^u, C M , C M , and & M gives rise to nontrivial Jacobians in the path integral, given by 



■h 



det' -□ 
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J c = [det'(-D)]- 1 , J b = [det'(-n)]- 1 



(18) 



These Jacobians can be absorbed by field redefinitions which however introduce terms which involve noninteger powers 
of the Laplacian. To avoid technical difficulties, we therefore prefer to exponentiate these Jacobians by the introduction 
of auxiliary anticommuting and commuting fields according to the sign of the exponent of the determinant, see also 
H, [l(j ■ One has to take their contribution into account while writing the FRGE. 



III. THE GAUSSIAN MATTER FIXED POINT 

The running of V a (<fi 2 ) is calculated from the FRGE as 

where (d t T k)[4> 2 , R] is obtained for various fields in an analogous way as in 0, [nj. Rescaling all fields with respect to 
the cutoff scale k, we obtain the dimensionless quantities <j) = k^4>, R = k~ 2 R and V a {4> 2 ) = k~( d ~ 2aS) V a {(j) 2 ). These 
dimcnsionlcss quantities we can use to analyze the RG flow and its FP structure. From the running of V a {4> 2 ) one 
can calculate the running of V a ((j> 2 ) using 

(d t Va)l4> 2 } = -(d-2a)V a (4> 2 ) + (d-2)4> 2 v:(4> 2 ) + k-^- 2 ^(d t V a )[^ 2 } (20) 

where the last term is calculated using eq. (fT^|) . A FP is a solution of the infinite set of functional equations dtV a = 
for a = 0, . . . , oo. This means that, at the FP, for each a the function V a (4> 2 ) is fc-independent, or equivalently that 
each coefficient of its Taylor expansion is fc-independent. Since we assume that each V a is analytic it can be Taylor 
expanded around (f> 2 = 0, and therefore 

9^«(0)=0 (21) 
for i = 0, . . . , co, where the superscript i denotes the i-th derivative with respect to <fi 2 . 



A. Minimal matter coupling of gravity at the GMFP 

The existence of a Gaussian Matter Fixed Point (GMFP), where all the matter couplings approach zero for k — > co 
and only the purely gravitational couplings have nontrivial values, was observed for finite polynomial truncations in 
Q . In [l| , its existence was proven for effective average actions of the form 



T k [g, <f>\ = J d d x Vff (V o (0 2 ) + V^cf 2 ) R+± d^d^ + S GF 



Igh- (22) 
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The existence of a GMFP can be shown to hold for the more general class of effective average actions considered in 
this paper. By definition, a GMFP is a point where V a are (j) 2 -independent, i.e. 

K«(0) = (23) 

for i = 1, ...,oo. In this subsection we will prove that with the ansatz in cq. (|2"3")l all the equations in (12"T]) with 
i = 1, . . . , oo are identically satisfied, thus leaving only the equations with i = to be solved. We will give numerical 
solutions of these remaining equations for a = 0, 1 . . . , 8 in section ITVl 

Now we explicitly analyze the structure of dtF related to the second variation of the effective average action given 
in eq. (J5]) for the various field components. The second variation for TnK v and ^ has the form 



r (2) 

1 k 



= f(z,R) + f a (z,R)V a , (24) 

T,V 



(2 

where we denote z := — □. The functional form for T k 



is motivated by eqs. (J9j) and (Ti"0|) from which we notice 

that it depends on V a at most linearly, with coefficients being functions of z and R, which are denoted here by f(z, R) 
and f a (z, R). 

For the scalar part, the second variation has the form 



F (2) 
1 k 



l 11 (z,R) + ti 1 (z,R)V a l 12 (z,R) + f^(z,R)V a g^z, R) <f>v; 
l 12 (z 7 R) + f 1 a 2 (z 7 R)V a l 22 (z,R) + f 22 (z 7 R)V a g 2 a (z,R)cj>V: (25) 

gl(z, R) <j> VI g 2 a (z, R) ^V' a z + R a (2 V' a + 4 4> 2 K") 



where a prime denotes derivative with respect to <p 2 . Again the functional form for is motivated by eq. (fTTj) 

s 

which clearly tells that entries S aa , S a h, and Shh depend at most linearly on V a , while the entries and S$h are 
linear combinations of 0V a '. The coefficients in these linear combinations are functions of z and R denoted here by 
l»(z,R) and gl(z,R). 

For the ghost part the second variation has the form 



r (2) 

1 k 



= D{z, R) . (26) 

gh 



This can be verified from eqs. (fT4l [T5l [T6l and [T7|) . We first consider the contributions from and Since for 

them the second variation has the form given by eq. (|24[) . the modified inverse propagator Vu '■= rj? + Ttk and the 
cutoff IZk will have the functional form 

V k = f(P k ,R) + f a (Pk,R)V a , TZ k = f(Pk,R) - f(z, R) + {f a (Pk,R) - f a (z, R)}V a , (27) 

where we have simply replaced z by P k (z) := z + Rk{z) to obtain the modified inverse propagator. R k {z) is a profile 
function which tends to k 2 for z — > and approaches zero rapidly for z > k 2 . The RG-timc derivative of the cutoff 
TZk in eq. (|27|) is 

d t n k = d t f(P k ,R) + dtf a (P k ,R) V a + {f a (P k ,R) - fa(z, R)} d t V a . (28) 
Using eq. (|28|) in the FRGE one finds that the contributions from h^ v and £ M have the form 

d t V a = H a {V c ) + H ab (V c )d t V b . (29) 



This can be justified by noticing that dtR-k given by eq. (|28|) depends at most linearly on dtV,. On the r.h.s. of the 
FRGE, dtR-k occurs in the numerator, while the denominator contains the modified inverse propagator given in eq. 
(|2"T|) which depends at most linearly on V a . So we find that the r.h.s of the FRGE depends at most linearly on dtV a . 
The coefficients in front of d t V a are functionals of V a and are denoted by H a (V c ) and H ao (V c ). 

The contributions from the ghost parts will be simpler. Since they do not depend on the potentials, they will only 
give a constant contribution to H a . The contributions from the scalars are more involved due to the matrix structure. 
The modified inverse scalar propagator is obtained by replacing all z with Pk in eq. 
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The cutoff is constructed in the usual way by subtracting the inverse propagator from the modified inverse propa- 
gator. This cutoff can be written as 

l n (P k ,R)-l u (z,R) l 12 (P k ,R)-l 12 (z,R) 
11% = ( l 12 (P k ,R)-l 12 (z,R) l 22 (P k ,R)-l 22 (z,R) 

P k - z 

flHPk,R)-fl 1 {z,R) f 1 a 2 (P k ,R)-f 1 a 2 (z,R) o\ 
/ 12 (P fc ,i?)-/ a 12 (z,i?) f 2 a 2 (P k ,R)-f 22 (z,R) V a 
/ 

gl(P k ,R)-gl(z,R) 

g 2 a (P k ,R)~g 2 a (z,R) | ■ (30) 

gl(P k ,R)-gl(z,R) g 2 a (P kl R) - g 2 a {z,R) 



Then the t derivative of the cutoff given in cq. (|30|) is 

d t l n (P k ,R) d t l 12 (P k ,R) \ / dtf^(P k ,R) d t fl 2 {P k ,R) 0' 
dtlll = | d t l 12 (P kl R) d t l 22 (P k ,R) + d t fl 2 {P k ,R) d t f 22 (P k ,R) | V a 
d t P k ) \ 

f^(P k ,R)-f^(z,R) f 1 a 2 (P k ,R)~f 1 a 2 (z,R) 0\ / d t gl(P k ,R) 

f 1 a 2 (P kl R)-f a - 2 (z,R)f 22 (P kl R)-f 22 (z,R)0)d t V a +\ d t g 2 a (P k ,R) | <^ Q ' 

/ V^OPfc^) d t g 2 a (P k ,R) 

^(P fc)J R)-^,iJ)\ 

g 2 a (P k ,R)-g 2 a (z,R) <$>d t V' a . (31) 

g 1 a (P k ,R)-g 1 a (z,R) g 2 a (P k , R) - g 2 a (z, R) / 

The modihed propagator for scalars is the matrix inverse of eq. (|25[) with z replaced by P k . It is given by 

(nr^j^ww ■ (32) 

where Adj ("P|l) denotes the adjoint of the matrix (V k ) (the matrix of cofactors). The determinant is a functional 
depending only on V a , 4> 2 V^V b , and 2V' a + 40 2 V^'. This can be easily derived from the modified inverse propagator 
obtained from eq. (f2"5")) . 

All entries of the adjoint of consist of cofactors, thus it has the form 

Adj {V s k ) = I A 21 (V a , <j> 2 VX, 2 K + 40 2 Vi!) A 22 (V a , <j> 2 V' a Vi 2 V' a + 40 2 V a ") A 23 (V a )<f> V' a , (33) 
^ 31 (K)<X A 32 (V a )<f>v: A™(V a ) j 

where each entry depends additionally on P k and R. In order to calculate the RG trace, we multiply {T , k )~ 1 with 
dt1Z k and then take the matrix trace. Doing this we note that </>V a ' is either multiplied with another </>V a ' or it is 
multiplied with <fidt V' a - So the scalar contribution to the FRGE has the form 

d t V a \ s = H a (y a , <f> 2 VX, 2 K + 40 2 Va) + H ab (V a ,<f) 2 V'X, 2 V' a + 40 2 V0) d t V b 

+H abc (V a , 4? V'X, 2 V'a + 40 2 va) <t? Vi d t V' c . (34) 

The contributions from the transverse traceless tensor and transverse vector can also be combined in the above 
expression to write the full FRGE contribution in the same way as above. Then dtF = R a dtV a - 

After having calculated the structural form for the running of V a (4> 2 ), we use it to calculate the dimcnsionlcss beta 
functional using eq. (|2T)|) . which gives 

(d t V a ){4> 2 } = -(d - 2a)V a + (d- 2)0 2 V> a ' + H a (k, 4> 2 K K, 2V' a + 40 2 V>") 

+H ab (v a , 4> 2 V' a VI, 2V' a + 40 2 va) {(d - 2b)V b - (d - 2)~4> 2 V b ' + (d t V b )[4> 2 ]) 

+H abc (V a , j, 2 K V b , 2V' a + 40 2 va) ~tfVl { (d - 2c) \>' - (d - 2) (>V C " + V^j + (d t V c )'{4> 2 } } .(35) 
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Inserting eq. (f2Tj) in eq. ((35)) we get the fixed point equation 

= -(d - 2a)V a + (d- 2)0 2 K' + H a (y a , 4> 2 V„ V b \ 2V' a + 40 2 V Q ") 
+H ab (V a , 4? V; V{, 2V' a + 4ft Vf) {{d - 2b)V b - (d - 2)^'} 

+H abc (V a , 4> 2 K V b , 2V' a + 40 2 V?) 4> 2 V b ' { (d - 2c) T>' - (d - 2) (> V c " + V?) } . (36) 

The above equation is identically satisfied when we take its Taylor expansion around <ft 2 = and use eq. (|23[) . For 
example, taking one derivative with respect to </> 2 gives 

= -(d- 2a)V> + (d - 2) U?V' a < + K'} + + /J" + 2 K'V d ' + 2 f rf "y c ') 

L J oV c "[P Vc V d) 

+ {(d- 2b)v b - (d - 2)0 2 ^'} ^ v;' + si ~ 2 H ^ r) (v: v' d + 2 v': v' d + ~e y d " v' c ) 

+4> 2 V b " {(d- 2c)V c - (d- 2)0 2 1/'} # afcc + 2 T> b ' {(d- 2c)f c ' - (d - 2)0 2 l7 c " - (d - 2)1/'} # afcc 
+0 2 H' {(d - 2c) V c (d 2)0 2 K'} f + ^ (Fj K' + ^ 2 Fj' F e ' + 2 f e " Fj) 

+ . ^ . (2Vj 2) + 4Vf > + 44> 2 VP)) . (37) 

Setting </> 2 = and using the GMFP conditions, we see that the right hand side will be zero. One can take successive 
derivatives to verify that this property indeed holds when higher derivatives are taken. The only equation which is 
not automatically solved in this way is the one where we evaluate eq. (f3"S| at 4> 2 = and use eq. (|23[) . This is just 
the FP equation for an f{R) theory with a single minimally coupled scalar. We will solve these equations in section 

im 



B. Linearized Flow around the GMFP 



The attractivity properties of a FP are determined by the signs of the critical exponents defined to be minus the 
eigenvalues of the linearized flow matrix, the so-called stability matrix, at the FP. The eigenvectors corresponding 
to negative eigenvalues (positive critical exponent) span the UV critical surface. At the Gaussian FP the critical 
exponents arc equal to the mass dimension of each coupling, so the relevant couplings are the ones that are power- 
counting rcnormalizable (or marginally renormalizablc) . In a pcrtubatively renormalizablc theory they arc usually 
finite in number. 

At the GMFP, the situation is more complicated as the eigenvalues being negative or positive do not correspond 
to couplings being relevant or irrelevant. In principle, at the GMFP the eigenvectors corresponding to negative 
eigenvalues get contributions from all the couplings present in the truncation, thus making it more difficult to find 
the fixed point action. Thus understanding the properties of the stability matrix around the GMFP becomes crucial. 

Therefore we now discuss the structure of the linearized flow around the GMFP. It is convenient to Taylor expand 
the potentials V a (cj) 2 ) as 

K(0 2 )-E A ^( fc )^' (38) 

where are the corresponding couplings with mass dimension d — 2a — i(d— 2). We are assuming a finite truncation 
with up to p powers of R, i.e. a going from to p, and q powers of 2 . In practice it has been possible to deal 
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with p < 8; as we shall sec, it is possible to understand the structure of the theory for any polynomial in 2 , so 
one could also let q — *■ oo. Rescaling these couplings with respect to the RG scale defines dimensionless couplings 



Xn"' = k d 2a l( - d 2 ^X^) and the corresponding beta functions /3 2 ' 



;(<*) 

The stability matrix is defined as 



(«) 



a) 



(> f K (l) (0) 



dX 



FP 



(39) 



FP 



Using the above definitions, numerical results tell that the stability matrix M has the form 

\ 



M 01 








M n 


M 12 








M 22 


M 23 








M 33 



V 



(40) 



where each entry is a (p + 1) x (p + 1) matrix of the form 



Mi. 



( ap$> 








d/3 { 2 p) 





(41) 



while p is the highest power of scalar curvature included in the action. It turns out that, 

M i:j = Vi > 1, Vj < i ; M„ = Vi, V j > (i + 1) . 
The various nonzero entries follow the same relations that were observed in In d dimensions they are 



where 



M ( 



on 



Mi, 



( SMr(O) r(0) 
A ' A 



8Mu P ) r(0) 
/ (5Mr(o) r(0) 



Mqi 



V 5M xf 



(d -2)il + Mqo 

■ ■ ■ SMr(0) r(p) \ 



l5-Mjb) ^(p) y 



(< + l)(2i + l)M i 



(42) 
(43) 



/-d0 
-(d-2) 

-(d-4) 

V 



(44) 



-(d-2p)J 



5Mr ip 



(45) 



Using the same arguments as in [l|, one can prove the above properties starting from cq. (|34[) neglecting dtV a and 
dtVa on the right hand side (corresponding to a one- loop approximation). Solving eq. (|34[) beyond that level would 
require solving a functional differential equation and would be beyond the scope of this paper. However, the results 
presented in the next section suggest that these relations should hold exactly. They are relations independent of the 
gauge choice, however the entries of Afoo and Mqi are gauge dependent. 

The physical nature of the relations among the eigenvalues can be understood from the difference between the 
GMFP and the Gaussian fixed point where also the gravitational couplings would vanish. At a Gaussian fixed point, 
the critical exponents are determined by the mass dimension of the couplings, and therefore arc all spaced by d — 2. 
At the GMFP, the gravitational couplings lead to some corrections to the critical exponents, but the correction to all 
exponents is the same, such that the spacing remains equal to d — 2. 
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These relations have important consequences. Because the stability matrix at the GMFP has the block diagonal 
structure given by eq. (jlU)) . its eigenvalues are just the eigenvalues of the diagonal blocks. Since the diagonal blocks are 

related by eqs. (|43|) . the eigenvalues of the various blocks differ only by multiples of d — 2. That means if p$ , . . . , p p ^ 
are the eigenvalues of Moo, then all the eigenvalues of M are of the form 

P® =P i a) + (d-2) l . (46) 

As Moo depends only on the couplings A Q a ' , it is enough to include only these couplings into the action to find all the 
eigenvalues of the stability matrix. Therefore, the results for minimally coupled scalar-tensor theory determine the 
eigenvalues of the nonminimally coupled scalar-tensor theory. In particular, if one has calculated the dimension of 
the UV critical surface of the minimally coupled theory, one can also predict the dimension of the UV critical surface 
of the nonminimally coupled theory. 

To find all the eigenvectors of the stability matrix it is necessary to know also M i. One can write the eigenvectors 
as v = (i>o, V\, . . . , v q ) T where each Vi is itself ap+1 dimensional vector. Then the vector Vq = (vq, 0, 0, ... , 0) T is an 
eigenvector if Vq is an eigenvector of Mqq which can be seen immediately by multiplying it with M. The eigenvectors 
of M with the above form are eigenvectors for the eigenvalues of Moo and can therefore be completely determined by 
just using Moo- Thus we note at this point that these eigenvectors are mixtures of gravitational couplings only, they 
do not contain any contribution from matter couplings. 

Now consider a vector of the form V\ = (v' , Vi, 0, 0, . . . , 0) T . Acting on it with M, and demanding V\ to be an 
eigenvector of M corresponding to some eigenvalue p^\ we obtain two relations, 

Moo v' + M i vi = P { 2 a) v' , Mil vi = P { 2 a) V! . (47) 

The second equation in (|47|) tells that v\ is an eigenvector of M\\. Now due to equations given in (|43|) and (j46]l . we 
note that V\ = Vq. Determining V\ will then determine also v' . In the same way one can go on to determine the next 
eigenvector. Consider V 2 = (v'q, v[, v 2 , 0, . . . , 0) T . We then demand it to be a eigenvector of M. That means it should 
satisfy 

Moo v'q + M i v[ = p^' v' ' , Mn v[ + Mi 2 v 2 = p^ v[ , M 22 v 2 = p^ v 2 ■ (48) 

One notices immediately that v 2 is the eigenvector of M 22 , and using equations in (|43f and (|46|) we conclude that 
v 2 = vq. Other equations would determine v' ' and v[. This process can be continued to find all the eigenvectors. 

We will now illustrate the validity of these results in various truncations with scalar fields coupled minimally and 
nonminimally to gravity. 



IV. NUMERICAL RESULTS 
A. Nonminimally coupled scalar field 

From here on we proceed as in [13] • We choose the gauge a = 0, (3 — > oo, and p = 0. This simplifies the calculation 
considerably because with that choice several arguments in the FRGE cancel with each other. The cutoff operators 
are chosen so that the modified inverse propagator is identical to the inverse propagator except for the replacement 
of z = —V 2 by Pk(z) = z + Rk(z); we use exclusively the optimized cutoff functions Rk(z) = (k 2 — z)9(k 2 — z) [27l ] . 
Then knowledge of the heat kernel coefficients which contain at most R 4 taken from [28[ is sufficient to calculate all 
the beta functions. A further benefit of this choice of cutoff is that the trace arguments will be polynomial in z. This 
simplifies the integrations in the trace evaluation and is done in closed form. 

Inserting everything into the FRGE and comparing the terms with equal powers of R and </> 2 on each side of the 
equation will give a system of algebraic equations for the beta functions of the couplings . The fixed points of the 

flow equations are evaluated and the corresponding critical exponents i? 2 i are determined. 

We carried out the calculation for effective average actions including up to R 4 and up to (f> 2 in each potential V a - 
Such truncations include at most ten couplings. We find that a GMFP does indeed exist for all these truncations. 

The nonvanishing fixed point values for various truncations are given in table [H the corresponding critical exponents 
(the negative of the eigenvalues of the stability matrix) in table [HI 

From the critical exponents one realizes at once several features. Though we carry out the full FRGE calculation 
we find that in general the real parts of the critical exponents z?2 differ from i?q exactly by two as proven in the 
one-loop case while the imaginary parts of the critical exponents are unchanged. This suggests strongly that the 
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p 


~ n 

A o* 


A o* 


A <)» 


r 3 
A o* 


A 0* 


1 


6.495 


-21.579 








2 


5.224 


-16.197 


1.834 






3 


6.454 


-20.756 


1.071 


-6.474 




4 


6.354 


-21.342 


0.792 


-6.807 


-3.865 



TABLE I: Nonvanishing couplings at the GMFP. The index p is the highest power of R included in the truncation. All values 
are multiplied by a factor 1000. 



p 


K 






v o 




#2 


#i 


*w 






1 


2.493 


2.368 








0.493 


2.368 








2 


1.826 


2.366 


21.822 






-0.174 


2.366 


19.822 






3 


3.077 


2.524 


2.033 


-3.852 




1.077 


2.524 


0.033 


-5.852 




4 


3.261 


2.772 


1.670 


-3.593 


-5.182 


1.261 


2.772 


-0.330 


-5.593 


-7.182 



TABLE II: Critical exponents at the GMFP. The index p is the highest power of R included in the truncation. Critical 
exponents are labeled v$ ■> u ^ e * ne couplings, but the corresponding eigenvectors involve strong mixing, as discussed in the 
text. For each i, the first two critical exponents form a complex conjugate pair given by i9q ± flgi and $ 2 ± $2*- 



relations among the eigenvalues will also hold at the exact level. The qualitative and quantitative properties turn out 
to be very similar to those of the purely gravitational theory. 

The inclusion of only four couplings with a = 0, 1 and i = 0, 1 leads to four attractive directions. The complex 
critical exponents $' Q ± d'^i are expected from the experience with the Einstein-Hilbert truncation. The existence of 
a second pair of complex critical exponents z? 2 ± ?? 2 i follows from the relation between the eigenvalues given in eq. 
(|46|) . These complex conjugate pairs occur also when higher scalar curvature terms are included. 

When one includes also R 2 couplings, one encounters large positive critical exponents as known from the calculations 
in pure gravity [H, [l(| HH, [23| . Using eq. (f4"6")) one concludes that one has to go up to power </> 20 before encountering 
a negative critical exponent, so the critical surface would be twelve dimensional. But this is a fluke of the R 2 
truncations due to the anomalously large positive critical exponent. The situation quickly normalizes when one adds 
further powers of R. 

Including i? 3 couplings, classically one would expect only three positive critical exponents as the classical mass 

0, —2, 2, 0, —2, and —4 respectively. Apparently, 



dimensions of An°\ \{ , A 2 X , A 2 2 \ and X\ 0> , are 4, 2, 











,(3) 



the FRGE calculation, which includes quantum corrections with large mixing between the various couplings, produces 
instead six positive critical exponents in the R 3 truncation. The critical exponent 1)^ is however very close to zero 
in consistency with the eigenvalue shift in eq. (|46p . This tells us that the truncation with p — 3 has a six-dimensional 
UV critical surface for any i > 1. 

With the inclusion of the coupling for the i? 4 operator whose classical mass dimension is —4, one notices that 
< i?[, 2) < 2. Thus one would expect that including the coupling for the operator <fi 2 R 4 with classical mass dimension 

—6, in consistency with eq. (|4"rj|) . the critical exponent ?9 2 would be negative, and the critical surface would be 

(2) 

five dimensional. Indeed, the inclusion of those couplings does make i9 2 negative, leading to five negative and five 
positive critical exponents. One can then say, using eq. (|46|) in the truncation p = 4, that for any i > 1, the critical 
surface would be five dimensional. 

To illustrate our results we display here the stability matrix for the R truncation. The entries in the upper left 
5x5 block and in the lower right 5x5 block are the same except the ones on the diagonals which differ by two. The 
upper right block is Mqx, the lower left one contains only zero entries: 



M\ 



GMFP 



-0.81 

-8.01 

2.16 

2.95 

5.12 













1.87 
-6.05 
0.27 
-0.61 
4.95 













0.40 
2.95 
-4.57 
-7.46 
3.34 













-1.24 
2.78 
1.64 
4.13 
-10.52 













0.41 
-1.80 
-0.041 
0.44 
7.79 













-0.0057 
-0.0031 
0.00021 
-0.00026 
0.00065 

1.19 
-8.01 

2.16 

2.95 

5.12 



0.0021 
-0.0093 
-0.00018 
-0.0032 
0.0021 
1.87 
-4.05 
0.27 
-0.61 
4.95 



0.0011 
0.00083 
-0.0032 
-0.0098 
-0.0010 
0.40 
2.95 
-2.57 
-7.46 
3.34 



-0.00039 
0.0024 
-0.00038 
-0.0019 
-0.0071 
-1.24 

2.78 

1.64 

6.13 
-10.52 



-0.000051 \ 
0.00024 
-5.5510-6 
-0.000091 
-0.00075 
0.41 
-1.80 
-0.041 
0.44 

9.79 / 



(49) 
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The eigenvectors corresponding to the five positive critical exponents in the i? 4 truncation are given by 



V 



-0.2774 ± 0.2693z\ 

0.8574 
-0.1206 ±0.0634z 
0.0473 ± 0.1254i 
-0.2202 ±0.1746i 





o / 



(15.381 ± 5.409i) x 10~ 4 
(-33.008 ± 13.931*) x 10~ 4 
(4.894 ± 1.980i) x 10" 4 
(-2.535 ±1.083i) x 10~ 4 
(5.437 ± 8.333i) x !0~ 4 
-0.2774 ± 0.2692z 

0.8574 
-0.1205 ±0.0634i 
0.0473 ± 0.1254i 
-0.2202 ±0.1746i 



V 



0.3845 \ 
0.07586 
-0.7103 
-0.5667 
-0.1437 





o J 



(50) 



The first complex conjugate pair of eigenvectors corresponds to the complex conjugate pair of critical exponents 
$0 ± i9q with values 3.2608 ± 2.7722i, while the second pair of complex conjugate eigenvectors corresponds to the 
complex conjugate pair of critical exponents 'd^i'd'^ with values 1.2608 ± 2. 7722z. The last eigenvector corresponds to 
the critical exponent # = 1.6698. We note that the eigenvectors corresponding to the eigenvalues of Moo, namely 
the first complex conjugate pair of eigenvectors and the last one, have the same structure as was described in the 
previous section, i.e. (vq, 0, 0, ... , 0) T , where vo is determined by just using Moo- We note that these eigenvectors do 
not get mixing from the matter couplings, but only from the purely gravitational couplings. Further more, if we look 
at the eigenvectors corresponding to the eigenvalues of Mn, namely the second complex conjugate pair of eigenvectors 
in eq. (|50p . which has the form (v' Q , vi, 0, . . . , 0) T , we clearly notice that v\ = vq, as described in the previous section. 



B. Minimally coupled scalar field 

Having verified that the properties of the stability matrix proved at one-loop level do also hold in the exact calcu- 
lation, we now analyze higher order curvature terms retaining only the couplings A °^ corresponding to a truncation 
with a minimally coupled scalar field. Then one obtains the non-Gaussian fixed points and critical exponents given 
in tables IIIII and HVl We analyze these results and use them to make predictions for the nonminimal truncation. 



V 


A* 


G, 


A*G* 


10 3 x 










A o» 


A o* 


XW 
A o* 


A o» 


XW 
A o. 


A o. 


A o* 


xm 

A o. 


\(») 
A o, 


1 


0.150 


0.923 


0.139 


6.495 


-21.579 
















2 


0.161 


1.228 


0.198 


5.224 


-16.197 


1.834 














3 


0.155 


0.958 


0.149 


6.454 


-20.756 


1.071 


-6.474 












4 


0.149 


0.932 


0.139 


6.354 


-21.342 


0.792 


-6.807 


-3.865 










5 


0.149 


0.932 


0.139 


6.355 


-21.339 


0.793 


-6.793 


-3.854 


-0.024 








6 


0.146 


0.918 


0.134 


6.312 


-21.669 


0.586 


-7.169 


-5.576 


-0.537 


2.702 






7 


0.146 


0.917 


0.133 


6.318 


-21.702 


0.534 


-6.469 


-5.530 


-1.979 


2.761 


2.565 




8 


0.148 


0.926 


0.137 


6.344 


-21.489 


0.678 


-5.922 


-4.574 


-2.074 


1.863 


2.393 


0.829 



TABLE III: Position of the FP for increasing number p of couplings included. The first three columns give the FP values in 
the form of cosmological and Newton constant and their dimensionless product. The values Aq"' (and only them) have been 
rescaled by a factor 1000. 

One observes that the addition of the scalar fields alters the results for pure gravity in @, HO] only by a small amount. 
Just as there, the UV critical surface becomes at most three-dimensional, and fixed point values for the cosmological 
and the Newton constant remain very stable. It has to be remarked that for those two couplings the oscillation in the 
fixed point value after the introduction of the i? 2 -term is not as strong as in pure gravity. Also the critical exponent 
obtained after the introduction of the i? 2 -coupling becomes large, but not as large as in pure gravity. So the addition 
of the scalar field seems to have already a little stabilizing effect on the i? 2 -truncation. The introduction of the R A and 
i? 5 -couplings leads to a second complex conjugate pair of critical exponents as soon as both couplings are included. 
Now it is easy to analyze how the dimension of the UV critical surface changes under the introduction of nonminimal 

matter couplings. In general, if a critical exponent is negative then v^i wm a l so De negative for all i > 0. From 
table [TV] we see that $ a) < for all a > 3, thus all ^ < for all a > 3 and i > 0. However, since 4 > $' > 2, 
using cq. (Titjl) we conclude that 2 > d' 2 > 0. This means that there are two more attractive directions. From table HVl 
one sees however that < i? 2 ' ) < 2 as soon as i? 4 is included, thus we do not obtain any other attractive directions. 
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p 


#'o 


K 


,4(2) 


,0(3) 


,Q( 4 ) 






i9 (7) 

^0 


7 q( g ) 
v o 


1 


2.493 


2.368 
















2 


1.847 


2.397 


21.031 














3 


3.077 


2.524 


2.033 


-3.852 












4 


3.261 


2.772 


1.670 


-3.593 


-5.182 










5 


2.777 


2.908 


1.795 


-4.176 


-4.196 


-6.764 








6 


2.841 


2.813 


1.386 


-4.000 


-3.798 


-5.947 


-8.538 






7 


2.930 


2.964 


1.312 


-4.009 


-2.760 


-4.623 


-7.459 


-11.166 




8 


2.331 


2.902 


1.570 


-4.063 


-0.673 


-7.120 


-7.323 


-9.854 


-11.611 



TABLE IV: Critical exponents for increasing number p of couplings included. The first two critical exponents are a complex 
conjugate pair of the form -ff' ± ■&" i. The same is the case for the fourth and fifth critical exponent $q ^ ^cf'*- 

So compared to [|| [To| where a three-dimensional UV critical surface was obtained for pure gravity, interactions with 
scalar matter lead to a five-dimensional UV critical surface. 

V. CONCLUSION 

We have shown that a Gaussian matter fixed point does exist also under the inclusion of higher order curvature 
terms and their coupling to scalar fields. We verified that the properties of the stability matrix proven only at one-loop 
level hold also in the exact calculations. We exploited these properties to show the relations between minimally and 
nonminimally coupled scalar-tensor theory. In particular, we were able to calculate the critical exponents for the 
nonminimal scalar tensor theory from those of the minimal one. The introduction of minimally coupled scalar matter 
fields gives only slight quantitative corrections to the fixed point properties of the purely gravitational theory. The 
critical exponents again seem to converge with the inclusion of more curvature terms. The minimally coupled theory 
produces three positive critical exponents. We derived that the additional critical exponents in the nonminimally 
coupled theory will be the ones of the minimal theory shifted by constant values. This produces two more positive 
critical exponents. From that we can conclude that, in four dimensions, the scalar-tensor theory based on an action 
polynomial in scalar curvature and in even powers of scalar field gives rise to a five-dimensional UV critical surface. 

Acknowledgements 

We would like to thank R. Percacci for many useful conversations and advice on the manuscript. C.R. would like 
to thank Daniel Litim for discussions, and Gutenberg-University, Mainz, and SISSA, Trieste for hospitality during 
different stages of this work. 



[1] G. Narain and R. Percacci, "Renormalization Group Flow in Scalar-Tensor Theories. I", arXiv:0911.0386>ep-th]. 

[2] S. Weinberg, Lectures presented at Int. School of Subnuclear Physics, Ettore Majorana, Erice, Sicily, Jul 23 - Aug 8, 1976. 

Published in Erice Subnucl.Phys. 1976:1; In General Relativity: An Einstein centenary survey, ed. S. W. Hawking and W. 

Israel, pp. 790-831; Cambridge University Press (1979). 
[3] S . Weinberg, larXiv:0903.0568l [hep-th] ; larXiv:0908.1964l [hep-th] . 
[4] G. 't Hooft and M. J. G. Veltman, Annales Poincare Phys. Theor. A 20 (1974) 69. 

[5] S. Deser and P. van Nieuwenhuizen, Phys. Rev. D 10 (1974) 411; Phys. Rev. D 10 (1974) 401; Lett. Nuovo Cim. 11S2 (1974) 
218 [Lett. Nuovo Cim. 11 (1974) 218]; Phys. Rev. Lett. 32 (1974) 245; S. Deser, H. S. Tsao and P. van Nieuwenhuizen, 
Phys. Lett. B 50 (1974) 491; Phys. Rev. D 10 (1974) 3337. 

[6] A. Codello, R. Percacci and C. Rahmede, Annals Phys. 324 (2009) 414. 

[7] R. Percacci and D. Perini, Phys. Rev. D67, 081503(R) (2003). 

[8] R. Percacci and D. Perini, Phys. Rev. D68, 044018 (2003). 

[9] L. Griguolo and R. Percacci, Phys. Rev. D 52, 5787 (1995) [arXiv:hep-th/9504092| . 
[10] A. Codello, R. Percacci and C. Rahmede, Int. J. Mod. Phys' A 23 (2008) 143 . 

[11] M. Reuter and F. Saueressig, JCA P 0509 (2005) 012 [arXiv:hep-th /0507167| ; M. Reuter and F. Saueress ig, JCAP 0509 
(2005) 012 [arXiv:hep-th/0507167| ; A. Bonanno and M. Reuter, JCAP 0708 (2007) 024 )arXiv:0706.0174l [hep-t h]]. 

[12] C. Cheung, P. Creminelli, A. L. Fitzpatrick, J. Kaplan and L. S enatore, JHEP 0803 (2008) 014 [arXiv:0709.0293l [hep-th]] : 
S. Weinberg, Phys. Rev. D 77 (2008) 123541 |arXiv:0804.429Tl [hep-th]]. 



14 



[13] F.L. Bezrukov, M. Shaposhnikov, Phys. Lett. B659 703-706 (2008) arXiv: arXiv:0710.3755 [hep-th]; A. De Simone, M.P. 

Hertzberg, F. Wilczek, Phys. Lett. B67 81-8 f2009) larX"iv:0812.4946l [hep-ph]: CP. Burgess, Hyun Mi n Lee, M. Trott, 

JHEP 0909:103 (2009) larXiv:0902.4465l [hep-ph]; T.E. Clark, Boyang Liu, S.T. Love, T. te r Veldhuis, larXiv:0906.5595l 

[hep-ph]; A.O. Barvinsky, A.Yu. Kamenshchik, C. Kiefer, A. A. Starobinsky, C.F. Steinwachs, arXiv:0910.1041 [hep-ph]. 
[14] C. Wett erich, Astron. A strophys. 301 321-328 (1995) e-Print: |hep-th/9408025| Phys. Rev. Lett. 102 141303 (2009), 

e-Print: larXiv:0806.0741 1 [hep-th], 
[15] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, JCAP 050 2 (2005) 010; S. Nojiri and S. D. Odintsov, 

Int. J. Geom. Meth. Mod. Phys. 4 (2007) 115; T. P. Sotiriou and V. Faraoni. la7Xiv:0805. 172 61 [gr-qc]. 
[16] I. Antoniadis, P.O. Ma zur and E. Motto la, Phys. Rev. Lett. 79 14-17 (199 7) e-Print: |astro-ph/9611208| Phys. Lett. B 444 

284-292 (1998) e-Print: |hep-th/ 9808070 New J. Ph ys. 9 11 (2007) e-Print: |gr-qc/0612068| E. Mottola and R. Vaulin, Phys. 

Rev. D 74 064004 (2006) e-Print: |gr-qc/0604051| 

[17] P. F. Machado and R. Percacci, Phys. Rev. D 80, 024020 (2009) [arXiv: 0904 .25101 [hep-th]]. 
[18] C. Wetterich, Phys. Lett. B 301, 90 (1993). 
[19] M. Reuter, Phys. Rev. D57, 971 (1998). 

[20] D. Dou and R. Percacci, Class. Quant. Grav. 15 (1998) 3449; W. Souma, Prog. Theor. Phys. 102, 181 (1999). O. Lauscher 
and M. Reuter, Int. J. Mod. Phys. A 17, 993 (2002); R. Percacci, Phys. Rev. D73, 041501(R) (2006); M. Niedermaier, 
JHEP 0212 (2002) 066; Nucl. Phys. B 673, 131-169 (2003); Phys. Rev. Lett. 103 (2009) 101303; D. F. Litim, AIP Conf. 
Proc. 841 (2006) 322; P. Fischer and D. F. Litim, AIP Conf. Proc. 861 (2006) 336; O. Lauscher and M. Reuter, Phys. 
Rev. D65, 025013 (2002); Class. Quant. Grav. 19, 483 (2002); M. Reuter and F. Saueressig, Phys. Rev. D65, 065016 
(2002); P. Fischer and D. F. Litim, Phys. Lett. B 638 (2006) 4 97; M. Reuter and H. Weyer, Phys. Rev. D 79 (2009) 105005 
; M. Reuter and H. Weyer, Phys. Rev. D 80 (2009) 025001 [arXiv:0804.1475l [hep-th]]; Gen. Rel. Grav. 41 (2009) 983; 
O. Zanusso, L. Zambelli, G. P. Vacca and R. Percacci, arXiv:0904.0938 [hep-th]. E. Manrique and M. Reuter, Phys. Rev. D 
79 (2009) 025008; arXiv:0905.4220 [h ep-th] ;laxXly :0907.2617] [gr-qc]; J. E. Daum, U. Harst and M. Reuter. larXiv:0910.4938l 
[hep-th]; J. E. Daum and M. Reuter. [arXiv:0910. 5401 [hep-th]; A. Eichhorn, H. Gies and M. M. Scherer. larXiv:0907.182"8l 
[hep-th]. 

[21] O. Lauscher and M. Reuter, Phys. Rev. D 66, 025026 (2002). 
[22] A. Codello and R. Percacci, Phys.Rev.Lett. 97, 221301 (2006). 

[23] P. F. Machado and F. Saueressig, Phys. Rev. D 77 (2008) 124045 . 

[24] D. Benedetti, P. F. Machado and F. Saueressig. larXiv:0901.2984l [hep-th]: larXiv:0902.4630l [hep-th], 

[25] Max Niedermaier and Martin Reuter, Living Rev. Relativity 9, (2006), 5; M. Niedermaier, Class. Quant. Grav. 24 (2007) 
R171; R. Percacci, In "Approaches to Quantum Gravity: Towards a New Understanding of Space, Time and Matter", ed. 
D. Oriti, Cambridge University. larXiv:0709.385T1 [hep-th] : D. F. Litim. larXTv:0810.3675l [hep-th] . 

[26] I.L. Buchbinder, S.D. Odintsov and I.L. Shapiro, "Effective action in quantum gravity", IOPP Publishing, Bristol (1992). 

[27] D. Litim, Phys. Lett. B 486 (2000) 92 ; D. F. Litim, Nucl. Phys. B 631 (2002) 128 ; Phys.Rev. D 64 105007 (2001), 
Phys.Rev.Lett. 92 201301 (2004). 

[28] I. G. Avramidi, Lect. Notes Phys. M64 (2000) 1. 



